This program examines tracts that only partially intersect with community districts.
library(sf)
## Linking to GEOS 3.9.0, GDAL 3.2.1, PROJ 7.2.1
library(tidyverse)
## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --
## v ggplot2 3.3.5 v purrr 0.3.4
## v tibble 3.1.4 v dplyr 1.0.7
## v tidyr 1.1.3 v stringr 1.4.0
## v readr 2.0.1 v forcats 0.5.1
## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()
library(magrittr)
##
## Attaching package: 'magrittr'
## The following object is masked from 'package:purrr':
##
## set_names
## The following object is masked from 'package:tidyr':
##
## extract
path_name <- "../data/"
load(paste0(path_name, "bl", ".RData"))
glimpse(bl)
## Rows: 39,903
## Columns: 4
## $ bl_id <chr> "200910518012010", "200910514003006", "200910515002005", "200~
## $ bl_name <chr> "Block 2010", "Block 3006", "Block 2005", "Block 1009", "Bloc~
## $ bl_area <dbl> 38500.823, 51711.677, 104537.047, 73933.762, 95854.528, 19934~
## $ geometry <POLYGON [°]> POLYGON ((-94.64213 38.9774..., POLYGON ((-94.63966 3~
load(paste0(path_name, "cd", ".RData"))
glimpse(cd)
## Rows: 59
## Columns: 4
## $ cd_id <dbl> 106, 108, 113, 102, 129, 116, 114, 101, 105, 103, 107, 109, 1~
## $ cd_name <chr> "East Side", "Old Northeast", "Greater Downtown", "Blue Valle~
## $ cd_area <dbl> 281713850, 118237219, 181920811, 216842678, 412671242, 264282~
## $ geometry <POLYGON [°]> POLYGON ((-94.52337 39.0941..., POLYGON ((-94.50777 3~
load(paste0(path_name, "icd", ".RData"))
glimpse(icd)
## Rows: 15,748
## Columns: 19
## $ bl_id <chr> "290950023002014", "290950063002017", "290950022003001",~
## $ bl_name <chr> "Block 2014", "Block 2017", "Block 3001", "Block 3004", ~
## $ bl_area <dbl> 14983.197, 9597.627, 18151.722, 18380.541, 17691.230, 17~
## $ st_id <chr> "29", "29", "29", "29", "29", "29", "29", "29", "29", "2~
## $ co_id <chr> "29095", "29095", "29095", "29095", "29095", "29095", "2~
## $ tr_id <chr> "29095002300", "29095006300", "29095002200", "2909500220~
## $ bg_id <chr> "290950023002", "290950063002", "290950022003", "2909500~
## $ people <chr> "62", "0", "22", "26", "42", "49", "43", "0", "12", "20"~
## $ housing_units <chr> "18", "0", "9", "18", "22", "22", "14", "0", "1", "10", ~
## $ land_area <chr> "14996", "9606", "18167", "18396", "17706", "17473", "32~
## $ water_area <chr> "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "~
## $ block_lat <chr> "+39.0894021", "+39.0439974", "+39.0932160", "+39.093298~
## $ block_lon <chr> "-094.5318307", "-094.5610343", "-094.5240738", "-094.52~
## $ cd_id <dbl> 106, 106, 106, 106, 106, 106, 106, 106, 106, 106, 106, 1~
## $ cd_name <chr> "East Side", "East Side", "East Side", "East Side", "Eas~
## $ cd_area <dbl> 281713850, 281713850, 281713850, 281713850, 281713850, 2~
## $ geometry <GEOMETRY [°]> POLYGON ((-94.5324 39.0887,..., POLYGON ((-94.5~
## $ in_area <dbl> 14983.19694, 9597.62653, 18151.72239, 18380.54058, 17691~
## $ pct_in <dbl> 100.0000000, 100.0000000, 100.0000000, 100.0000000, 100.~
load(paste0(path_name, "red", ".RData"))
glimpse(red)
## Rows: 39,903
## Columns: 11
## $ st_id <chr> "20", "20", "20", "20", "20", "20", "20", "20", "20", "2~
## $ co_id <chr> "20091", "20091", "20091", "20091", "20091", "20091", "2~
## $ tr_id <chr> "20091050000", "20091050000", "20091050000", "2009105000~
## $ bg_id <chr> "200910500001", "200910500001", "200910500001", "2009105~
## $ bl_id <chr> "200910500001000", "200910500001001", "200910500001002",~
## $ people <chr> "52", "47", "51", "58", "138", "68", "35", "81", "39", "~
## $ housing_units <chr> "33", "29", "27", "26", "70", "34", "18", "28", "12", "1~
## $ land_area <chr> "27826", "26904", "26385", "27075", "74489", "58378", "2~
## $ water_area <chr> "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "~
## $ block_lat <chr> "+39.0427053", "+39.0427081", "+39.0426993", "+39.042717~
## $ block_lon <chr> "-094.6217603", "-094.6228101", "-094.6238297", "-094.62~
load(paste0(path_name, "bg", ".RData"))
glimpse(bg)
## Rows: 1,684
## Columns: 4
## $ bg_id <chr> "201039819003", "201030714004", "201030711041", "201030712051~
## $ bg_name <chr> "Block Group 3", "Block Group 4", "Block Group 1", "Block Gro~
## $ bg_area <dbl> 984975.5, 36602869.0, 18313082.5, 26007918.3, 11448660.8, 243~
## $ geometry <POLYGON [°]> POLYGON ((-94.94756 39.3373..., POLYGON ((-95.18756 3~
icd %>%
tibble %>%
distinct(cd_id) %>%
arrange(cd_id) %>%
pull(cd_id) -> cd_list
for (i in cd_list) {
icd %>%
filter(cd_id==i) %>%
tibble %>%
distinct(bg_id) %>%
arrange(bg_id) %>%
pull(bg_id) -> bg_list
print(
paste(
i,
paste(
bg_list,
collapse=", "),
sep=": "))
cd %>%
filter(cd_id==i) -> cd_i
bg %>%
filter(bg_id %in% bg_list) -> bg_i
for (j in bg_list) {
bg_i %>% filter(bg_id==j) -> bg_j
ggplot(bg_i, aes()) +
geom_sf(
fill=NA,
color="lightgreen") +
geom_sf(
data=cd_i,
aes(),
fill="white",
color="darkred") +
geom_sf(
data=bg_j,
aes(),
fill="NA",
color="darkgreen") -> map1
plot(map1)
}
}
## [1] "101: 290950078022, 290950105001, 290950105002, 290950105003, 290950105004, 290950105005, 290950106001, 290950106002, 290950106003, 290950106004, 290950107021, 290950107022, 290950107023, 290950119003, 290950119004, 290950120001, 290950120002, 290950120003, 290950124003, 290950124004, 290950125011, 290950125012, 290950125021, 290950126001, 290950126002, 290950128032, 290950128043, 290950129033, 290950156001, 290950163001, 290950163002, 290950163003, 290950170001, 290950170002, 290950171004, 290950172001, 290950172002, 290950172003, 290950172004, 290950174001, 290950174002, 290950174003, 290950175001, 290950175002, 290959801011, 290959808021"
## [1] "102: 290950019001, 290950019002, 290950019003, 290950020001, 290950020002, 290950021001, 290950021002, 290950021003, 290950021004, 290950021005, 290950022001, 290950022002, 290950022003, 290950023001, 290950034001, 290950034002, 290950058011, 290950058013, 290950110013, 290950155001, 290950155002, 290950155003, 290950156001, 290950156002, 290950156003, 290950163001, 290950163003, 290950164001, 290950170001, 290950171001, 290950171002, 290950171004, 290959808021"
## [1] "103: 290950072001, 290950073011, 290950073012, 290950074001, 290950074002, 290950074003, 290950074004, 290950075002, 290950075003, 290950081002, 290950082001, 290950082002, 290950082003, 290950082004, 290950083001, 290950083002, 290950086001, 290950086002, 290950086003, 290950086004, 290950086005, 290950086006, 290950087002, 290950090003, 290950091001, 290950091002, 290950091003"
## [1] "105: 290950061001, 290950063002, 290950075001, 290950075002, 290950075003, 290950076001, 290950076002, 290950076003, 290950077001, 290950077002, 290950078021, 290950078022, 290950079001, 290950079002, 290950079003, 290950080001, 290950080002, 290950080003, 290950081001, 290950081002, 290950087001, 290950087002, 290950088001, 290950171002, 290950171003, 290950171004, 290950175002, 290959801011"
## [1] "106: 290950022001, 290950022002, 290950022003, 290950023001, 290950023002, 290950034001, 290950034002, 290950034003, 290950034004, 290950037001, 290950038001, 290950038002, 290950052001, 290950053001, 290950053002, 290950054001, 290950055001, 290950056011, 290950056012, 290950056021, 290950056022, 290950057001, 290950057002, 290950058011, 290950058012, 290950058013, 290950060001, 290950061001, 290950061002, 290950061003, 290950063001, 290950063002, 290950075002, 290950077001, 290950154013, 290950158001, 290950160001, 290950160002, 290950161001, 290950161002, 290950162001, 290950162002, 290950164001, 290950165001, 290950165002, 290950166001, 290950166002, 290950169001, 290950169002, 290950169003, 290950171001, 290950171002, 290950171003"
## [1] "107: 202090452004, 290950046001, 290950046002, 290950046003, 290950051001, 290950051002, 290950052001, 290950053002, 290950063001, 290950065001, 290950065002, 290950066001, 290950066002, 290950067001, 290950067002, 290950069001, 290950069002, 290950071001, 290950071002, 290950071003, 290950071004, 290950072001, 290950072002, 290950073011, 290950073012, 290950073021, 290950073022, 290950074001, 290950074002, 290950075002, 290950167001, 290950167002, 290950167003, 290950167004, 290950168011, 290950168012, 290950168013, 290950168021, 290950168022, 290950169001, 290950178001, 290950178002, 290950178003, 290950178004"
## [1] "108: 290950003001, 290950003003, 290950006001, 290950006002, 290950006003, 290950006004, 290950007001, 290950007002, 290950007003, 290950008001, 290950008002, 290950008003, 290950008004, 290950009001, 290950009002, 290950009003, 290950009004, 290950010001, 290950010002, 290950010003, 290950018001, 290950018002, 290950018003, 290950018004, 290950019001, 290950019002, 290950019003, 290950019004, 290950023001, 290950154011, 290950154012, 290950154013, 290950154021, 290950154022, 290950155001, 290950160001, 290950161002"
## [1] "109: 290950087001, 290950087002, 290950088001, 290950088002, 290950088003, 290950088004, 290950089001, 290950089002, 290950090001, 290950090002, 290950090003, 290950090004, 290950090005, 290950094002, 290950095001, 290950095002, 290950095003, 290950095004, 290950096001, 290950096002, 290950096003, 290959801011"
## [1] "110: 290950078022, 290950079001, 290950088004, 290950089001, 290950096002, 290950096003, 290950163002, 290950163003, 290950172002, 290950172003, 290950174003, 290950175001, 290950175002, 290959801011"
## [1] "111: 200910509004, 200910515004, 290950090003, 290950090004, 290950091001, 290950091002, 290950091003, 290950091004, 290950092001, 290950092002, 290950092003, 290950093001, 290950093002, 290950093003, 290950094001, 290950094002, 290950094003, 290950094004, 290950095002, 290950095003, 290950095004, 290950097002, 290950098001, 290950098002, 290950099001"
## [1] "112: 200910508001, 200910509001, 200910509002, 200910509004, 290950072001, 290950072002, 290950073011, 290950073021, 290950073022, 290950074003, 290950074004, 290950083001, 290950084001, 290950084002, 290950084003, 290950085001, 290950085002, 290950085003, 290950085004, 290950086003, 290950086004, 290950086005, 290950086006, 290950091003, 290950092001, 290950092003"
## [1] "113: 202090430004, 202099805001, 290470221006, 290950003001, 290950003002, 290950003003, 290950011001, 290950011002, 290950043001, 290950043002, 290950044001, 290950046001, 290950046002, 290950052001, 290950152001, 290950152002, 290950152003, 290950153001, 290950153002, 290950153003, 290950154011, 290950154012, 290950155001, 290950157011, 290950157012, 290950157021, 290950157022, 290950157023, 290950158001, 290950158002, 290950159001, 290950159002, 290950161001, 290950161002, 290950162001, 290950167003, 290950178001, 290950178004"
## [1] "114: 202099800001, 290470202012, 290470202021, 290470202022, 290470221002, 290470221006, 290950003002, 290950152002, 290950152003, 290950155001"
## [1] "116: 290470206022, 290470206024, 290470221006, 290470222003, 290950003002, 290950003003, 290950006001, 290950007002, 290950008003, 290950009001, 290950010003, 290950110011, 290950110012, 290950110013, 290950155001, 290950155002, 290950156002, 290950156003"
## [1] "117: 290470222003, 290470223013, 290470223023, 290950151001, 290950155002"
## [1] "118: 290470205001, 290470205002, 290470205003, 290470205004, 290470205005, 290470205006, 290470205007, 290470206021, 290470206022, 290470206023, 290470206024, 290470206025, 290470206031, 290470206032, 290470206033, 290470206034, 290470206041, 290470206042, 290470206043, 290470206044, 290470206045, 290470221004, 290470222001, 290470222002, 290470222003, 290470223023, 290950155001"
## [1] "119: 290470202024, 290470203001, 290470203002, 290470203003, 290470203004, 290470204001, 290470204002, 290470204003, 290470205001, 290470205007, 290470206033, 290470206034, 290470206044, 290470209011, 290470209012, 290470209013, 290470209014, 290470209015, 290470209021, 290470209022, 290470210014, 290470210032, 290470210033, 290470210044, 290470211044, 290470212103, 290470213072"
## [1] "120: 290470202011, 290470202012, 290470202013, 290470202014, 290470202015, 290470202021, 290470202022, 290470202023, 290470202024, 290470203004, 290470203005, 290470204002, 290470204003, 290470205001, 290470205005, 290470221001, 290470221002, 290470221004, 290470221006"
## [1] "121: 290470202022, 290470203005, 290470210011, 290470210012, 290470211011, 290470211015, 290470211016, 290470212041, 290470212042, 290470212043, 290470212044, 290470212081, 290470212082, 290470212083, 290470212084, 290470212085, 290470212086, 290470212141, 291650300022, 291650300023, 291650300031, 291650300041, 291650300042, 291650300043, 291650300044, 291650301011, 291650301012, 291650301021, 291650301022, 291650301023, 291650301031, 291650301033, 291650301034, 291650302012, 291650302013, 291650302071, 291650302072, 291650302073, 291650302074, 291650302075, 291650302076, 291650302081, 291650302082, 291650302101, 291650302111, 291650302112, 291650302113, 291650302121, 291650302122, 291650302123, 291650302131, 291650302132, 291650302141, 291650302151, 291650302161, 291650302162, 291650303081, 291650303082, 291650303083, 291650303084"
## [1] "122: 290470208021, 290470208022, 290470208031, 290470209011, 290470209021, 290470210033, 290470211034, 290470211044, 290470212101, 290470212102, 290470212103, 290470212104, 290470212121, 290470212123, 290470213031, 290470213032, 290470213033, 290470213071, 290470213072, 290470213073, 290470213074, 290470213075, 290470213076, 290470213091, 290470213092, 290470213093, 290470213111, 290470213112, 290470213121, 290470213122, 290470213123, 290470213143, 290470218071, 290470218072, 290470218081, 290470218082, 290470218083, 290470218102, 290470219003, 290470219004"
## [1] "123: 290470211011, 290470211012, 290470211031, 290470211032, 290470211033, 290470211034, 290470212081, 290470212083, 290470212085, 290470212091, 290470212092, 290470212093, 290470212101, 290470212102, 290470212111, 290470212112, 290470212121, 290470212122, 290470212123, 290470212131, 290470212132, 290470212141, 290470213031, 290470213032, 290470213033, 290470213131, 290470213132, 290470213133, 290470213134, 290470213141, 290470213142, 290470213143, 290470213144, 290470219003, 290470219004, 291650302013, 291650302132, 291650306012"
## [1] "124: 290470212141, 290470219003, 291650302011, 291650302012, 291650302013, 291650302014, 291650302082, 291650302101, 291650302132, 291650303062, 291650303081, 291650303083, 291650303084, 291650303101, 291650306011, 291650306012, 291650306025, 291659800001"
## [1] "125: 200910516002, 200910517002, 200910517003, 290950093002, 290950095002, 290950096002, 290950096003, 290950097001, 290950097002, 290950098001, 290950098002, 290950098003, 290950098004, 290950099001, 290950099002, 290950100011, 290950100012, 290950100021, 290950100022, 290950100023, 290950101031, 290950101032, 290950101033, 290950102011, 290950102012, 290950172003, 290950173004"
## [1] "126: 290950101032, 290950102011, 290950102012, 290950102031, 290950172003, 290950174003"
## [1] "127: 290950102011, 290950102031, 290950129031, 290950129032, 290950129033, 290950129034, 290950129041, 290950129042, 290950129043, 290950129062, 290950129063, 290950129065, 290950130031, 290950130032, 290950130033, 290950130034, 290950131003, 290950144001, 290950174001, 290950174002, 290950174003"
## [1] "128: 290950128022, 290950128023, 290950129042, 290950129065, 290950136061, 290950143002, 290950144001, 290950144002, 290950176004"
## [1] "129: 290950124002, 290950124003, 290950125021, 290950125023, 290950127021, 290950127031, 290950127032, 290950128022, 290950142031, 290950142033, 290950142051, 290950142062, 290950143001, 290950143002, 290950143003, 290950143004, 290950144002, 290950145011, 290950145012, 290950145013, 290950145041, 290950145042, 290950176001, 290950176002, 290950176003, 290950176004"
## [1] "130: 200910517003, 200910533011, 200910533013, 200910533023, 290950100012, 290950100022, 290950101031, 290950101032, 290950101033, 290950101051, 290950101052, 290950101053, 290950101054, 290950102033, 290950102044, 290950134011, 290950134012, 290950134052, 290950134071, 290950134072, 290950134171, 290950134172, 290950134173, 290950134181, 290950134182, 290950134183, 290950173001, 290950173002, 290950173003, 290950173004"
## [1] "131: 290950101032, 290950102011, 290950102031, 290950102032, 290950102033, 290950102034, 290950102041, 290950102042, 290950102043, 290950102044, 290950130032, 290950131003, 290950133014, 290950134011, 290950134183, 290950174003"
## [1] "132: 290950102031, 290950102042, 290950129061, 290950129063, 290950129064, 290950129065, 290950130031, 290950130032, 290950130033, 290950131001, 290950131002, 290950131003, 290950132031, 290950132032, 290950132033, 290950132034, 290950132081, 290950132082, 290950132101, 290950132102, 290950132103, 290950133014, 290950133071, 290950133072, 290950144001, 290959883001"
## [1] "133: 290370600044, 290370603072, 290370603073, 290370603091, 290370603092, 290950129065, 290950133131, 290950134161, 290950135043, 290950136061, 290950136124, 290950136141, 290950144001, 290950181021, 290950182001, 290959883001"
## [1] "134: 200910533024, 200910534092, 200910534093, 200910534291, 200910534311, 200910534312, 290370600011, 290370600013, 290370600041, 290370600042, 290370601001, 290370601002, 290370601003, 290950134012, 290950134051, 290950134052, 290950134072, 290950134101, 290950134103, 290950134161"
## [1] "201: 202090402001, 202090409001, 202090411001, 202090419002, 202090429005, 202099800001, 202099805001, 202099809001, 290470221006, 291650301012"
## [1] "202: 202099800001, 202099805001, 202099809001, 290470221006, 290950152003"
## [1] "203: 202090401002, 202090401003, 202090402001, 202090406003, 202090407001, 202090407002, 202090409001, 202090411001, 202090412001, 202090412002, 202090413002, 202090414001, 202090414002, 202090414003, 202090416002, 202090429001, 202090429004, 202090429005, 202090445002, 202099800001"
## [1] "204: 202090411001, 202090419001, 202090419002, 202090420011, 202090429001, 202090429003, 202090429004, 202090429005, 202099800001"
## [1] "205: 202090401001, 202090401002, 202090401003, 202090405001, 202090405002, 202090406001, 202090406002, 202090406003, 202090407001, 202090412001, 202090413001, 202090413002, 202090445002"
## [1] "206: 202090416002, 202090416003, 202090419001, 202090419002, 202090420011, 202090420012, 202090420021, 202090420022, 202090421001, 202090421002, 202090421003, 202090422001, 202090422002, 202090423001, 202090423002, 202090423003, 202090424001, 202090424002, 202090424003, 202090429001, 202090429002, 202090429003, 202090429005, 202099809001, 202099812001"
## [1] "207: 202090406001, 202090406003, 202090412001, 202090413001, 202090413002, 202090413003, 202090413004, 202090415001, 202090415002, 202090416001, 202090416002, 202090416003, 202090416004, 202090422001, 202090422002, 202090423003, 202090429002, 202090439032, 202090439052, 202099812001"
## [1] "208: 202090415002, 202090419002, 202090424001, 202090424003, 202090426001, 202090426002, 202090426003, 202090430001, 202090439052, 202099805001, 202099809001, 202099812001, 202099815001"
## [1] "209: 200910500001, 200910500002, 200910500003, 200910501001, 200910501002, 202090430002, 202090430003, 202090430004, 202090433011, 202090433012, 202090433013, 202090451001, 202090451002, 202090451003, 202090452001, 202090452002, 202090452003, 202090452004, 202090452005, 290950046002, 290950071001, 290950071002, 290950153002, 290950168021, 290950168022"
## [1] "210: 202090428001, 202090428002, 202090430001, 202090430002, 202090430003, 202090430004, 202090433011, 202090433012, 202090451002, 202090452004, 202099805001, 202099809001, 202099812001, 202099815001, 290950152003, 290950153002"
## [1] "211: 200910501001, 200910502004, 200910503013, 200910522011, 202090427001, 202090427002, 202090428001, 202090428002, 202090428003, 202090430001, 202090430002, 202090433011, 202090434001, 202090434002, 202090435001, 202090435002, 202090436004, 202090437001, 202090437002, 202099815001"
## [1] "212: 202090427001, 202090428001, 202090437001, 202090438022, 202090438031, 202090439052, 202090440024, 202099812001, 202099815001"
## [1] "213: 200910522011, 200910522012, 200910522013, 200910522021, 202090436001, 202090436002, 202090436003, 202090436004, 202090436005, 202090437001, 202090437002, 202090438022, 202090438023, 202090438031, 202090438032, 202090438033, 202099815001"
## [1] "214: 200910522022, 200910523051, 200910523061, 200910523062, 202090436001, 202090438021, 202090438022, 202090438023, 202090438031, 202090440023, 202090440024, 202090447023, 202099815001"
## [1] "215: 202090439051, 202090440011, 202090440012, 202090440013, 202090440021, 202090440022, 202090440023, 202090440024, 202090447021, 202090447023, 202090447042"
## [1] "216: 202090438022, 202090439052, 202090440013, 202090440021, 202090440022, 202090440023, 202090440024, 202090447023"
## [1] "217: 202090405002, 202090413001, 202090413004, 202090415001, 202090415002, 202090439031, 202090439032, 202090439041, 202090439042, 202090439051, 202090439052, 202090440013, 202090440024, 202090441011, 202090441041, 202090444003, 202099812001, 202099815001"
## [1] "218: 202090439041, 202090439051, 202090440011, 202090440012, 202090440013, 202090441011, 202090441021, 202090441022, 202090441031, 202090441032, 202090441033, 202090441041, 202090441042, 202090441043, 202090442013, 202090442021, 202090447021, 202090447031, 202090447042"
## [1] "219: 202090401003, 202090402001, 202090405001, 202090443013, 202090443021, 202090444001, 202090444002, 202090445001, 202090445002, 202090446012, 202090446031"
## [1] "220: 202090442021, 202090442022, 202090443011, 202090443012, 202090443013, 202090443022, 202090446011, 202090446012, 202090446014, 202090446021, 202090446031"
## [1] "221: 202090402001, 202090445001, 202090445002, 202090446011, 202090446012, 202090446021, 202090446031, 202099800001, 291650301012, 291650301032, 291650303064"
## [1] "222: 202090405001, 202090405002, 202090406001, 202090413001, 202090439031, 202090439032, 202090439042, 202090443021, 202090443031, 202090444001, 202090444002, 202090444003, 202090445002"
## [1] "223: 202090439041, 202090439042, 202090442022, 202090443011, 202090443012, 202090443013, 202090443021, 202090443022, 202090443031, 202090443032, 202090446012"
## [1] "224: 202090439041, 202090441031, 202090441032, 202090441041, 202090441042, 202090442011, 202090442012, 202090442013, 202090442014, 202090442021, 202090442022, 202090442023, 202090443022, 202090443032, 202090447041, 202090447042"
## [1] "225: 202090442011, 202090442012, 202090442021, 202090446011, 202090446013, 202090446014, 202090446021, 202090447031, 202090447032, 202090447041, 202090447042, 202090448041, 202090448042, 202090448043"
## [1] "226: 201030711023, 201030711041, 201030712021, 201030712042, 201030712051, 202090446021, 202090447021, 202090447031, 202090447032, 202090447033, 202090448031, 202090448032, 202090448033, 202090448034, 202090448041, 202090448042, 202090448043"
## [1] "227: 201030711023, 202090446011, 202090446014, 202090446021, 202090446031, 202090448041, 202090448043, 291650303064, 291650304014"